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Abstract 

We obtain an explicit simple formula for the coefficients of the 
asymptotic expansion for the factorial of a natural number, 



n>. nV2^e-"{l + ^ + ^+^ + ---}, 

n n 



in terms of derivatives of powers of an elementary function that we call 
normalized left truncated exponential function. The unique explicit ex- 
pression for the afc that appears to be known is that of Comtet in [6] , 
which is given in terms of sums of associated Stirling numbers of the 
first kind. By considering the bivariate generating function of the as- 
sociated Stirling numbers of the second kind, another expression for 
the coefficients in terms of them follows also from our analysis. Com- 
parison with Comtet's expression yields an identity which is somehow 
unexpected if considering the combinatorial meaning of the terms. It 
suggests by analogy another possible formula for the coefficients, in 
terms of a normalized left truncated logarithm, that in fact proves to 
be true. The resulting coefficients, as well as the first ones are identi- 
fied via the Lagrange inversion formula as the odd coefficients of the 
inverse of a pair of formal series. This in particular leads to the iden- 
tification of a couple of simple implicit equations, which permits us to 
obtain also some recurrences related to the ajJ s. 

Keywords: Y function, asymptotic expansions, Lagrange inversion formula, Stirling 
numbers. 
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1 Introduction 



Consider, for t > 

which means that we may think e~* j x as the probability pj that a random 
variable X t takes the value j. This amounts to say that X t has a Poisson 
distribution with parameter t, whose expectation, variance and characteristic 
function, (p x (9) for 9 G ffi are: 

E(X t ) = t Var(X t ) = t 

^(9) := E(e»*) = Y,*"'* = E e "^ = ^^"^ ( L1 ) 

Take now t = n G N, and consider the random variable 

X n - n 



n 



;i. 2 ) 



From the Central Limit Theorem we know that the characteristic function 
of Z n converges pointwise to that of a standard Normal random variable: for 
9eR, 

E{j dZn ) ^e~ e2/2 asrwoo. (1.3) 

On the other hand, projecting the series (II. ip with t = n over e~ ien in 
L 2 [— 7r,7r] we obtain 



/TV /* TV — n 71 

d 0E(e» x «)e-» n = / d9 Q -^- = 2vr- 
-7T J -IT n - 



TV. 



Recalling the definition ( II. 2p of Z n , a change of variables in the first integral 
above yields 



d9E(e iez ") = / d9E{e l ^ {Xn - n) ) = 2n 



SO 



n\ = V ^ e " nnn (1.4) 
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If the convergence in (jl.3p could be seen to hold uniformly, we would recover 
Stirling's formula: 

n! « V27me~ n n n 

It is not difficult to see that this is indeed the has been noticed by 

several authors, see for instance [12] and the references therein, or Sect 27, 
exercise 18 of [5]. L. Baez-Duarte considered in [2] the above method to 
obtain the asymptotics for the general term of a convergent power series of 
a function / with non-negative terms, under conditions that guarantee that 
a uniform Central Limit Theorem for the corresponding random variables 
holds and, as a particular instance, he obtains the asymptotics for p(n), the 
number of partitions of an integer n. The approach has a longer history, and 
goes back to the articles by W.K. Hayman [9 J and E.A. Bender [3J. 

Here, we obtain an asymptotic expansion for the denominator in (jl.4p 
in terms of inverse powers of n, and show that it can be explicitly inverted, 
yielding a simple expression for the coefficients at of the expansion of n!, in 
terms of an elementary function (see Theorem 12.11 below) . To our knowledge, 
the unique explicit formula for those coefficients is that given by Comtet 
as an exercise in page 267 of [6], in terms of sums of derangements d% of 
integers with certain restrictions, related with the Stirling numbers of the 
first kind. As a corollary of the proof of our first formula, after recalling 
the bivariate generating function for the associated Stirling number of the 
second kind S3, we obtain a second expression for the coefficients, which is 
exactly alike to that in [5J, and in particular yields an unexpected equality 
between certain alternating sums of associated Stirling numbers of the first 
and second kind. This identity suggested another expression for the a&, in 
terms of a normalized left truncated logarithm function corresponding exactly 
with the first one, that we proved to be true (see Theorem 12. 7p . Finally, we 
recognize in these two expressions the odd terms of the series obtained by the 
Lagrange inversion formula applied to appropriate functions. Each of these 
two inverses satisfy a simple implicit equation. By differentiation we obtain 
a couple of differential equations that lead to a pair of recurrent relations 
that could be used also to generate the a^'s. One of them is similar to that 
obtained in [TT] using other techniques (see also [7]). 

Although we state the results as expressions for the coefficients of the 
expansion of n\ in inverse powers of n, they could be stated as well as expres- 
sions for the coefficients of the asymptotic expansion of the gamma function 
in inverse powers of z G C, for \z\ — > 00 in the region arg(z) G [0, 7r). This 
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follows at once from the existence of this last expansion. For this and other 
properties of the gamma, see PQ. Let us also note that the first few coeffi- 
cients, that in particular appear in most books and tables of special functions 
(for instance [8 J or [10j), can be computed one by one by exponentiating the 
well known Stirling series for the logarithm of the gamma, or from recurrence 
relations obtained by different techniques (as in [TJ] or [H]). We think how- 
ever that having explicit and simple formulae has an intrinsic interest. In 
particular, we show that they render some combinatorial identities, and re- 
late with the Lagrange inversion formula for elementary functions, what helps 
to underestand the relation of the expansion with the Lambert W function 
appearing in [7]. 

We provide precise statements concerning the asymptotic expansions, and 
obtain the formulae for the coefficients in the next section. In the last sec- 
tion, we show the relationship with the Lagrange inversion formula, obtain 
recursive formulae and generalise the identities obtained between sums of S3 
and d 3 . 



2 The asymptotic expansion and formulae 
for the coefficients 

Let us denote by d k f the kth derivative of a function /, with respect to its 
real or complex variable. Our first result is the following formula for the 
coefficients of the expansion of n! in powers of -. 

Theorem 2.1. Let n6N. Then the coefficients of the expansion 



n! 

fc>0 

are given by a k = -L 9 2fc (G-^)(0) (2.2) 

e x — 1 — x \ ~\ x- 7 
where G[x) = 2 = 2 > — 

Remark 2.2. Formula (12.11) is to be understood as an asymptotic expansion, 
that is, for any given N > 0, 

n\ = \P2^ e~ n n n ( ^ — a k + R N+ ^j 

k=0 
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with Rn+i = Q( n iv+i ) as n — > oo. See for instance Jlty or [L$ on the subject 
of asymptotic expansions. 

The function G above is what we call a normalized left truncated expo- 
nential. In general we define 

Definition 2.3. Given a function F (x) = fn^r with f 2 7^ 0, its normalized 
left truncated associated function F 2 is defined as 

F(x) -h-hx 

We also obtain another expression for at in terms of the 3— associated 
Stirling numbers of second kind. Recall that, for r an integer > 1, n,k 
integers > 0, the r— associated Stirling number of the second kind S r (n,k) 
is the number of partitions of a set of n elements into k blocks, all with at 
least r elements. The convention is that SV(0,0) = 1. Then we have 

Theorem 2.4. The coefficients in (12.11) are also given in terms of the 
^-associated Stirling numbers of second kind by 

VV 1V - S 3 (2(j + k),j) 
ak = 2^(j + k)\ (2 ' 3) 



Proof of Theorem \2.1\ The strategy will consist in expanding the denomina- 
tor in ( II. 4p in powers of -, and then taking the inverse of the resulting series. 
From flU]) and (JO}, for 9 G R, 

E(e i6Z ") = e^ 6 *- 1 "^) = e ^- e2 ^\ (2.4) 

where, for z G C 



e z - 1 - z - z 2 /2 



z 2 ^(k + 2)V 

fc>i v ' 



For given K G N, consider a Taylor expansion of order K in powers of z of 

c" 



K 



e -■»(*) 



J2jW(e-° 2 ^)\^ + R K+1 (z) (2.6) 



j=0 
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The remainder Rk+x satisfies the estimate 

I ~\K+l 



\Rk+i(z)\ < 



— sup \d K+1 (e- e2 ^)\ 



(2.7) 



Substitution of (12. 4p and (12. 6p in the denominator D n in (II. 4p yields 



D, 



2vT J —iTy/n 



d9E(e WZn ) 



2vT J —iTy/n 



d6e 2 e yy ^' 



'2ix 



j=0 



2 = 



+ 



due 2 



(2.8) 



V27T J-nyfE V 71 

It can be easily seen by induction that, for j G N and z G C 

where Pj, y (-) is a polynomial of degree j whose coefficients depend on prod- 
ucts of derivatives of g of order up to j evaluated in z. From ( 12. 51) it is easy 
to see that g and all its derivatives are bounded over compacts. Moreover, if 
we denote by TZ(z) the real part of z, we have 



(2.9) 



and, for x G M, 

K(g(ix)) =K 



e lx — 1 — ix + x 2 /2 



— ar 



cosx — 1 + x 2 /2 



x' 



X 



— — d 4 cos (ax) for some a G [0,1] 



(2.10) 



From (12. 7ft . (12. 9 ft and (12. 10ft . we can estimate the last integral in ( 12.81) as 
follows: 

d9e~~ R K (—=] 



< 



d6e~^ 



§1 | id \K+1 f)2nl i //l2 \l - Ck 



n 1 



e y ^sup|Px + i,,(^)|<^r, (2.11) 



■ <7T 



n 2 
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where Ck is a constant which depends on K but not on n. Then, observe 
that the difference of each integral in the sum in (12. 8p with respect to the 
integral in the whole R is exponentially small in n: for any j > 0, there is 
some a positive, 



x=0 



< 



d0e" 



|6»|>ttV^ 



x=0 
2 



< e" 



(2.12) 



Moreover, if we interchange the derivative with the integral in the second 
line above, and then compute the integral, we obtain: 



1 



d9e~ 



a/tv ?! V / 



x=0 



(^I(i-l)!!^((l + 2^))-^) 



x=0 



if j odd 
if j even 



(2.13) 



Then, for any given positive N G N, consider K = 2N + 1 in (12. 8p . sum and 
substract the integrals in the whole line to each term in the sum, estimate 
the differences using (I2.12p . evaluate the integrals as in (I2.13P and estimate 
the term with the remainder using (12. lip . Rename finally the summation 
index to obtain: 

N (—D k 1 1 
foranyiV>0, D n = £ L_L __ ^( G -^)( ) + 0(-_), (2.14) 

where we have used that 1 + 2g(x) = G(x), as defined in ( 12. 2p . Let us 
shorthand (I2.14p above as 



and recall we have 



fc>0 



(2.15) 



(2.16) 
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Consider next the the well known Stirling series, which is an asymptotic 
expansion for large values of \z\ of the logarithm of the gamma function in 
powers of \ (see for instance Theorem 1.4.2 in [I]). Its coefficients are given 
in terms of the Bernoulli numbers Bk- In particular, it implies 

log( n !) x log (^e-„») + £ 2t(at 
Since the series S n above contains only odd powers of -, it follows that 



y2wne n n n 



log( j ) x -S n = S- n . 

Therefore, 

n k>0 

and the theorem follows from (12.1 6p □ 



To prove Theorem 12.41 we will use the generating function of the r- 
associated Stirling numbers already introduced. Namely, (see page 222, ex- 
ercise 7, Ch. 5 of [6]) : for any r > 1 

H(t,u) := e " ( ^ + (^ + - ) = S r (l,k)u kt ^ (2.17) 

l,k>0 



Proof of Theorem 2J^_ Observe that, from the first equation in (12. 41) we 
may write the integral in the denominator in (II. 4ft as 



d9E(e iez ") = / d0e-^e nf( %\ (2.18) 

2 

with f(z) = e z — 1 — z — y- If we expand 



e nf(z) = Y 4^e n/ (0), 
i>o J ' 
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substitute in (I2.18P with z = and proceed as in the proof of Theorem 12. 1^ 
(that is, consider a finite Taylor expansion with remainder, interchange the 
sum and derivative with the integral, and observe that the integral may be 
considered in the whole line with an error smaller than the remainder) we 
obtain that for any given K, 

K_ -j 1 poo 2 

The remainder can be seen to satisfy Rr+i = 0{ (JC + 1)/2 ). After computing 
the integral and renaming the terms, we may shorthand the above as 

Next, from the generating function for S3 (I2.17P and observing that Ss(N, k) = 
if k > N > 0, we obtain 

^V'(0) = 5> fe S 3 (2j,A;) 

k<2j 

Substituting in (12.191) . and then summing over I = k — j, after recalling that 
from the combinatorial interpretation of S3 it is easy to see that S 3 (2j, k) = 
if A; > j > 

V ( v (-i)^ (-i) fc ^(2(fc + /),fc) 

^ > n j~k 2 jj\ 2^ n i 2^ 2 k + l (k + i)\ 

j>0,k<2j J l>0 k=0 v ' 

Concluding as in the proof of 12. II that = -D_ n , we obtain (12. 3 p □ 



From Theorems 12.11 and 12.41 since the asymptotic expansion of a function 
is unique, we have that 

1 ^/..Mw . ^ ( ■ S 3 (2(j + k),j) 
7Fk\ d {G ) (0) = 2^(j + k)\ 

Moreover, we have a second identity that follows from the expression for 
the coefficients of the expansion for the 7 function given in [6], and already 



9 



mentioned. It is 



v fc>o i=o u ' 

Permutations without fixed points (in other words, without cycles of length 
1) are known as derangements, and for natural r > 1, d r (n, I) is the number 
of derangements of a set of n elements, that have I cycles, all of length > r. 
Equating the coefficients in (12.31) and the corresponding ones resulting from 
( I2.20p when considering n\ = nT(n), we get the remarkable combinatorial 
identity presented in the following proposition. 

Proposition 2.5. For S3 and d% as defined above and for any given positive 
k G N, 

* S 3 (2(j + k),j) » d 3 (2(j + fc),j) 

2 j+k U + ky. f^ 2i+ k (j + k)\' K ' 

This identity is somehow surprising, since it is needed a precise balance 
to control the different growth of S3 and d%. We give a direct proof of a 
generalization of it in Proposition 13.21 Now, let us compute the generating 
function for ds(j, k). 

Lemma 2.6. The generating function for d% as defined above is 

$(t,u) = (l-t)-« e -«C+£) = ^(hk)u k - v (2.22) 

Proof We first compute the exponential generating function of the cyclic 
permutations of length at least 3 



T,( n - i y--=J2-= i °^ i - t y 1 -( t +^ 

z — ' n\ z — ' n \ 2 



k>3 fc>3 

By the exponential formula (see [13], Corollary 5.1.6), we obtain that 

= (i_t)-i e -(* + ^) 
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is the exponential generating function for the number of derangements of a 
set of n elements, whose cycles are all of order at least 3. Then, the generating 
function 

$(t, u) = <j)(t) u = (1 - t)- u e ~ u{t+t ^\ 

also keeps track of the number of cycles on this kind of derangements, (see 
[T3] . Example 5.2.2) obtaining equation f[2T22j) . □ 

By writing the generating function H of S3 (see (I2.17P ) and $ in the form 
H(t,u) = A* 1 - 1 -^) = e «(-^-'M-4) 

the analogy between them is apparent, and, together with the identity (I2.2ip . 
suggests that a formula for the a k in terms of a normalized left truncated 
logarithm, in the sense of definition 12.31 may hold. This proves to be true, 
yielding another formula, exactly alike to (12. ip . 

Theorem 2.7. The coefficients in the asymptotic expansion of n\ 

n\ = a/27™ e~ n n n ^ ^ — ^ a k 
k>o n 

are also given by a k = d 2fc (L~^)(0) (2.23) 

2 fc k\ 

t t( \ o ~ + x ) + x r xj 

for L(x) = 2 - = 2 > — — - 

Proof Consider the identity 

POO 

T(x) = x x dte^H*- 1 x > 



After integrating by parts once, and changing variables y/x (t — 1) — > u, we 
obtain 

= x a e-V2^ /l r dte -*{t-i) e *iozt^\ (224) 

^e" 1 ^/ 1 Z" 00 , -x(-s-) xiog(-*.+i) 

awe V -A e V s 



e 
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Define the function 



to obtain, using the expansion 



■j 



in f l2T24j) with y 



i x = ; / due 2 e ^ 



27T J-^/i 

1= f° rfwe-^V (^V'IdV^(O) (2.26) 



j>0 

It is clear that 



\y=0 



2+1 



f j_ r ^^e-^C 1 -^)^ « H 1 - 2 ^))~ 2 ^' - if j even 

Vv 7 ^ i-v^ ' |0 oth 



2tt J-y^E ' I otherwise, 

(in the sense that the difference is exponentially small in i as i -)• oo, the 
details are analogous to that in the proof of [2. 1[) . Renaming and reorganizing 
terms in the sum above, and observing that 1 — 2£(y) = L(y) as defined in 
(T22SD it follows 

fc>0 

Taking x = n and using that n! = nT(n), we conclude the proof, from ( 12. 241) 
and the above formula for /. □ 

Since a proof of the expansion f 1 2 . 2 j) of Comtet does not seem to be 
available, we provide next a brief sketch of it, that follows from a slight 
modification of the arguments in the proof above, exactly as Theorem 12.41 
was obtained by modifying the proof of Theorem 12.11 
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Proof of HPS} . Write the integral in (EMD as 

J?/, e 2 +x where h(n) = lof fl -4- tA — n 4- - 

and expand 



I(x) = . / due "a" where = log (1 + y) — y + — . 

V27T y. ,/x 2 



e x/.&o = ^|_ 5 v fe (o). 



fc>0 

Substituting in the integral with y — 



i r 00 2 7/ fc 



k>0 

2 



fc>0 

and integrating as before 



dwe 2 u 



(A; — 1)!! if k even 



y27T 10 otherwise, 

after renaming and collecting terms, we have 

U x ) = Y —^r-d 2j e xh (0) 

Now, from ( 12. 221) and the definition of /i, <&(y,x) = e~ xh ^~ y \ and 

^e^(O) = d 2 'e xh ^\ y=0 = dPe xh (-\ =0 = d? j ${-x,y)\ y=0 
Hence, from (ET221) 

i>o J ' j>o J ' fc=o 



A; 



which changing k + I = j and proceeding exactly as in the last part of the 
proof of Theorem 12.41 yields ( I2.20P . 
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3 Lagrange inversion and recursive formulae. 



We will identify the ak's in (12.21) and (I2.23[) as factors of the odd coefficients of 
a pair of formal power series, with the aid of the Lagrange inversion formula, 
that we recall next in a suitable form. Let R(x) be a formal power series 
of the form R(x) = xT(x), where T(x) is a formal power series with non 
zero constant term. Denote by R^ 1 ^ the substitutional inverse of R, i.e., 
S = R^-V means S(x)T(S(x)) = S(xT(x)) = x. Let S(x) = J2T=i ^Jf- The 
Lagrange inversion formula gives a simple recipe to compute the coefficients 
of the series S (see for example jl]). 

Sk = d^ 1 (T~ 1 ) k (0) = d^ 1 (T- k ) (0). (3.1) 
Define the following exponential formal power series 

x , x k ( ,e x — 1 — Xs i/2\ 

B M=5>*r : =(*(-P72-) ) < 3 ' 2 > 

k=l 1 

= {(2e x - 2 - 2x) 1/2 ) { - 1) 

/t=i ■ ' 

= ((-21og(l + x) + 2x) 1/2 ) < - 1) . 

Observe that from equations (12. 2 p and ( 12.23jl . and the Lagrange inversion 
formula (13. ip we identify 

where bk '■— bk/k\ and Ck '■— c^/kl. Then, the coefficients of the formal power 
series B{x) and C(x) coincide at odd powers. A much stronger result holds 
indeed. 

Proposition 3.1. The formal power series B(x) andC(x) differ only at the 
coefficient of the quadratic term. More precisely, we have 

C(x) - B(x) = 



14 



Proof. From equations ( 13. 2 p and (13. 3p we obtain the implicit equations 

e B(x) _ 1 _ B ^ = ^ ^ 3 4 ) 



and 



x 2 



C( x )-- = \ g(l + C(x)). (3.5) 



Taking log at both sides of (13 .4p , 

2 

B(x) = \og (l + B(x) + ?j), (3.6) 

2 

which is equivalent to ( 13. 5 p by the change C{x) = B{x) + All the opera- 
tions to obtain respectively (13. 6 p and ( 13. 5 p from (13. 2p and (13. 3p are reversible 
in the context of formal power series. Then, the result follows. □ 

Explicit expressions for the coefficients of the series B(x) and C(x) can 
be obtained by using Lagrange inversion formula, leading to the following 
generalization of the identity (12.211) . 

Proposition 3.2. The following identity holds for every k ^ 2 

k-l 



S 3 (k + 2j-l,j) 



(* + l)(* + 3)...(* + 2j-l) 



^fc+j-i rf 3 (fc + 2j-l,j) 

^ l J (* + l)(* + 3)...(* + 2j-l)" 



Proof. By the Lagrange inversion formula we have 



b k = d 



k-l 



1 — X 



x 2 /2 



-k/2 



2 = 

( 1 + r- 1 -^^ 1 ^ ~ m 
X 2 



x=0 
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By the binomial identity 



^ + 2 e*-l-x-x 2 /2 ^ " fe/2 _ g ^/^j 2 i (e" - 1 - x - x 2 /2)^ 

oo 

= ^(-l)^ + 2)...(fc + 2(j-l)) /!2; 



(e x - 1 - x - x 2 /2)^' 
^— x;- a^/v -r ±) ■ ■ ■ \r>, t *>yj — *-)) 

3=0 



= + 2)...(k + 2(7 - !))^ E 5 s(U)|- (3.7) 

j=0 X Z>3j 

The last equality follows from equation (12.171) for r = 3, after writing the 
exponential in power series of the exponent, and equating powers of it, to 
obtain 

J ' l>3j 



Interchanging sums in (13.71) we get 

jk + 2)... 

/(/-!).. .(/-2j + l) ^ J> {l-2j)Y 



■ fc(fc + 2)...(fc + 2(j-l)) 
Z^Z^i- 1 ) f/_9ui) ^'^Tfr^r (3 " 9) 



Z=0 3j<l 

The coefficient of x k ~ x /{k— 1)! of this series, obtained by making l—2j = k—1, 
is equal to 

1V - g 3 (fc + 2j-l,j) 
2^ J (* + l)(* + 3)...(* + 2j-l) 



Notice that from equation 12.221 proceeding as to obtain (13. 8 ft 

(-log(l -x) - a; - x 2 /2) j v^,/,.^' 
_^d 3 (l,j)j. 

J ' l>3j 

From that we easily obtain 

(log(1+l) - x+x72)J - = E( _ i) ,^ (i , f . (3io) 
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By expanding the binomial 

log(l + x) -x + x 2 /2\ k/2 _ /-log(l +x) +x x k/2 
x~ 2 J ~\ ^72 

using equation (13.101) . and proceeding exactly as for the computation of bk, 
we get 

-V( nfc+i-i d 3 (k + 2j-l,j) 
k jrC ) (k + i)(k + 3)...(k + 2j-iy 

The result then follows from Proposition 13.11 □ 

Now we state and prove recursive formulae to obtain b^, b^, cj., and c^. 

Proposition 3.3. The sequences b k , c^, bk and c^ satisfy the following re- 
cursive formulae 



A- 



j 

k-2 



k = - jttt ( ~ hk - x + E (j + ~ hk - j ) (3 - 12) 

i=i 



Ck 



Y^n y k Ck - 1 - 12 Q) c i+i c ^j ( 3 - 13 ) 



fc-2 



( a 7— — ( Cfc-i - ^(j + 1) c j+ i4-i ) (3.14) 

J'=l 



fc + 1 



wift &! = Ci = &i = Ci = 1. 

Proof. Formulae (13.121) and (13.141) follow from formulae (13.111) and (13.131) 
respectively. Computing the derivative in both sides of equations ( 13.60 and 
( 13.50 . after clearing up we obtain 

B'(x)B(x) = x-{x 2 /2)B'{x) (3.15) 
C'{x)C{x) = xC{x)+x. (3.16) 

Recurrence (13.110 (respectively ( 13.130 ) is obtained by equating coefficients of 
x k /k\ in both sides of the resulting series in ( 13.151) (respectively ( 13.160 ). □ 
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Observe that ( I3.14p is similar to the recurrence obtained in [TTJ (see also 
[7]). We list a few terms of both formal power series B(x) and C(x), 

B(x) = x-x 2 /6 + a; 3 /36 -x7270 + x 5 /4320 + x 6 /17010 + ... 
C(x) = x + a; 2 /3 + x 3 /36-a; 4 /270 + x 5 /4320 + x 6 /17010 + .... 

The coefficients a^'s can be readily computed as well from formulae (I2.2p or 
(I2.23p . using for instance MAPLE. In particular, we checked that the first 
twenty terms coincide with that of [15] . 
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